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Abstract 

We continue our systematic construction of Baxter Q-operators for spin chains, which is based on 
certain degenerate solutions of the Yang-Baxter equation. Here we generalize our approach from 
the fundamental representation of gl(n) to generic finite-dimensional representations in quantum 
space. The results equally apply to non-compact representations of highest or lowest weight 
type. We furthermore fill an apparent gap in the literature, and provide the nearest-neighbor 
Hamiltonians of the spin chains in question for all cases where the gl(ra) representations are 
described by rectangular Young diagrams, as well as for their infinite-dimensional generalizations. 
They take the form of digamma functions depending on operator- valued shifted weights. 



1 Introduction, Overview and Outlook 



We continue our systematic approach to the construction of Baxter Q-operators for those types of 
quantum spin chains that are based on rational R-matrices |T] , [2] , [3] . Our program is motivated 
by the desire to gain a deeper understanding of the integrable structure of the free/planar 
AdS/CFT system, where spin chains appear in the weak coupling limit. We will, however, not 
directly explore this connection in this paper, and refer the reader to the remarks made in [I], 
[3] as well as to the overview article series |1] for further information. Instead, we will proceed 
to extend the methodology of [2], where only the fundamental representation was considered, 
to the case where the state space of the quantum spin chain is allowed to be a tensor product 
of arbitrary finite-dimensional or non-compact Ql(n) representations of highest or lowest weight 
type. Note that non-compact lowest-weight representations are precisely the ones which are 
needed in the AdS/CFT spectral problem, albeit in a supersymmetric generalization. The latter 
will be the subject of future work. 

In this article we restrict ourselves to the purely "bosonic" Ql(n) case, and provide in Section[2] 
a very general construction of the Baxter Q-operator. In Section [2.11 we briefly recall our earlier 
work, where we had algebraically found a new class of Lax operators for the fundamental repre- 
sentation of gl(n) by studying certain singular solutions of the Yang-Baxter equation. We also 
set up the generalization to general representations in quantum space. Section 12.21 proceeds to 
find the needed, more general solutions of the Yang-Baxter equation, yielding novel (quantum) 
Lax operators for arbitrary gt(n) representations. Of special importance in our construction 
is the non-commutative Cayley-Hamilton theorem, which is reviewed in Appendix lAl We then 
succinctly discuss the fusion properties of the Lax operators in Section 12.31 Up to this point 
our work is purely algebraic, and not restricted to specific representations. After discussing in 
12.41 certain analytic properties of the derived new class of Lax operators, we proceed, as in our 
earlier work, in Section [231 to construct Baxter operators by first building suitable monodromies 
based on the new solutions, and subsequently tracing over the oscillator representations of the 
auxiliary space, a technique pioneered in |5]. At this point we choose to restrict ourselves to 
lowest or highest weight representations, deferring the most general case to future work. As 
shown in Section I2.6| functional relations follow in the usual way from the fusion relations, and 
systems of Bethe equations are then deduced in a straightforward manner from these relations. 
We thereby rederive the most general system of nested gl(n) Bethe equations, first obtained in 
[6 J by means of the Bethe ansatz (which our method entirely avoids). A closely related method 
to obtain these equations called "analytical Bethe ansatz" was presented in [TJ. See also [8], 
where this generic system of Bethe equations was also derived (including the supersymmetric 
case) from the fusion relations of the eigenvalues of Baxter operators, without providing, how- 
ever, an operatorial construction. Our construction has the additional advantage that whenever 
Q-operators are diagonalizable then the set of solutions of the Bethe equations is complete. To 
conclude our Q-operator construction, we show in Section 12.71 and in particular in Appendix iBl 
that our method is also of very practical value, if one e.g. is interested in explicitly constructing 
Q-operators for chains of small length. This is done by providing Mathematica™ programming 
code. 

Section [3j which may be read largely independently, is concerned with providing, tongue-in- 
cheek, interesting problems to the solutions found in Section [2l Baxter operators and transfer 
matrices for integrable spin chains form a large family of commuting operators. For applications, 
one is often interested in a special member of this family: a local, nearest-neighbor Hamiltonian. 
Integrable Hamiltonians for general q [(2) representations have been known for a long time, [9], 
|10| . |11| . Their derivation is based on the general method of finding a suitable transfer matrix 
with a special "regular point" in the spectral parameter plane, where the matrix becomes the shift 
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operator on the chain [12J. The Lax operator (R-matrix) used in the construction of these special 
transfer matrices intertwines two identical representations of the type present in the quantum 
space of the spin chain. For gl(2) it takes the form of a ratio of Euler gamma functions, whose 
arguments depend on the joint quadratic Casimir operator of the intertwined representations. 
The Hamiltonian density of the spin chain is then obtained as the logarithmic derivative of 
the R-matrix w.r.t. the spectral parameter at its shift point value, and turns into (a sum of) 
digamma functions. See in particular |13j for a beautiful discussion of this approach. To our 
knowledge, an explicit generalization of these g[(2) formulas for the Lax operator and associated 
nearest-neighbor Hamiltonian density to the case of general gl(n) spin chains is lacking in the 
literature. The state of the art appears to be the so-called "tensor product graph method" |14| . 
|15| (and further references therein) for multiplicity-free representations, which expresses the 
R-matrices as a (finite or infinite) series involving projectors on all irreducible representations 
occurring in the two-site tensor product decomposition. We should also mention two other 
approaches to the construction of R-matrices. For finite-dimensional representations of fjl(n), 
one can also obtain R-matrices by applying the standard fusion procedure to Yang's R-matrix, 
see e. g. [16J. Furthermore, for continuous series representations of fll(rt) R-matrices are given 
in a remarkable factorized form in [17] . The earlier constructions [B], [7] do not propose the 
gl(re) generalization of the explicit Ql(2) Hamiltonians. This prompted us to (partially) fill this 
apparent gap in the literature, and to derive Hamiltonians for all those §l(n) highest or lowest 
weight representations that satisfy a "generalized rectangularity condition" (this includes the 
cases where the two-site tensor product is multiplicity- free) . We set up the derivation in Section 
13.11 and then solve in Section l3~2l the Yang-Baxter equation for R-matrices intertwining two such 
rectangular representations. Technically, the method very closely follows the one employed in 
Section [2 The result is expressed as products of gamma functions, whose arguments are linear 
expressions involving the spectral parameter and operator valued shifted weights. This explicit 
result is checked against the formal expression known from the "tensor product graph method" 
in Appendix O In Section \3. 31 we find the shift point of these R-matrices, where they turn into 
the permutation operator. Taking a logarithmic derivative, we find the Hamiltonian densities 
of quantum spaces in arbitrary rectangular representations as digamma functions functionally 
depending on the shifted weights. Some details are delegated to Appendix ID1 In Section \3. 41 we 
find, using the functional relations, the dispersion laws allowing to extract the spectrum of the 
Hamiltonians from the solution to the system of Bethe equations derived in Section 12.61 This 
is in line with our above statement that we provide a natural diagonalization problem (local 
Hamiltonians and their dispersion laws) to the solution (the systems of Bethe equations)! 

Let us mention some other recent work on the Baxter operators of fll(n) spin chains. An 
alternative method employing the so-called co-derivative was proposed in |18j . It would be 
interesting to see whether it may be generalized from the fundamental representation to general 
highest/lowest weight representations as in our current work. Furthermore, there are the articles 
[19J , [2D] ( see also many references therein for the general method), where the case of compact 
and non-compact £|[(2) representations is treated by an alternative procedure. It would be 
interesting to see whether this approach allows for an equally unified description of compact and 
lowest /highest non-compact gl(n) representations, as is needed when applying the Q-operator 
construction to e.g. the AdS/CFT system. 

We have obtained all of our new Lax operators for Q-operators (Section [2]) and for the 
transfer matrices with regularity property for spin chains in general representations (Section [3]) 
by explicitly solving the YBE on a case-by-case basis. It would be very interesting and, perhaps, 
conceptually more satisfying to recover all these results from Drinfel'd's universal R-matrix, see 
e. g. [21j, [22j for some related, recent work in this direction. This, as well as the generalization 
to the supersymmetric case, is left to future investigation. 
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2 Q-operators 



2.1 Review and Outset 

The quantum mechanical state space V of an integrable, homogeneous lattice chain model usually 
takes the form 

V = V ®V ® ...®V (2.1) 

* v ' 

L-timcs 

of an L-fold tensor product of equal "local" representation spaces V of some algebra. Here L is 
the length of the chain, and V is frequently abbreviated to "quantum space". In this article we 
take the algebra acting on V to be the Lie algebra fll(n), and will, as is customary, frequently 
refer to the lattice model as a "spin" chain, even though this name strictly speaking only applies 
to the case n = 2. 

Historically, integrable spin chains were discovered through solving, by means of Bethe's fa- 
mous ansatz, the diagonalization problem of some given, specific, nearest-neighbor Hamiltonians 
acting on V. The quantum inverse scattering method (QISM), see |13| for a very authoritative 
review, instead put the emphasis on the rather systematic construction of a commuting family of 
transfer matrix operators, which may be constructed by starting from the Yang-Baxter equation 
(YBE), and depend on a free diagonalization-independent parameter, the spectral parameter. 
The integrable nearest-neighbor Hamiltonian is then derived from this systematic approach by 
expanding the "correct" transfer matrix operator around the "correct" value of the spectral pa- 
rameter. See also Section [3] for a more detailed discussion. 

Curiously, while the usual transfer matrix operators forming the commuting family depend 
on the spectral parameter and an additional label, the representation in so-called auxiliary space 
|13| . the family is even bigger, and also encompasses certain singular transfer matrices termed 
Q-operators. Here "singular" carries a double sense: for one, according to the usual QISM these 
operators diverge unless the global 0l(n) symmetry of V, which for periodic spin chain boundary 
conditions is inherited from the g[(n)-invariance of V, is broken by certain regulating "angles". 
For another, these operators are constructed from certain singular solutions of the YBE. All 
this is discussed in our earlier work [2], where we constructed Q-operators for the integrable 
spin chain with V being in the fundamental representation of fll(n). In the following we extend 
our analysis to the general case where V is the representation space of any Ql(n) irreducible 
representation. 

According to the idea presented in [2J we construct Q-operators as special transfer matrices, 
where the auxiliary space is taken to be an appropriately chosen Yangian algebra representa- 
tion. Before proceeding to the determination of the relevant R-operator we briefly review the 
representations of the Yangian algebra introduced in |2J. The YBE 

R(*i - z 2 )(L{z 1 ) <8> 1)(1 <8> L(jzjj)) = (1 <g> L(2jj))(L(«i) <g> l)R(zi - z 2 ) (2.2) 

with 

R(z) = zl + P (2.3) 
provides the defining relation for the Yangian algebra y = Y(gl(n)). The L-operators of the 
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(2.5) 



satisfy the relation (|2.2p with L(z) = L/(z) provided that J h a form a subalgebra and 

(a^,a|) are families of oscillators, namely 
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d I — °dJb ~ °b Jd > 
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(2.6) 
(2.7) 
(2.8) 



Here / is any subset of the index set {1, . . . ,n}, and / is its complementary set. For graphical 
convenience the L-operator in (|2.4p is presented in the case where the set / = 
Any other L/(z) can be obtained by appropriate permutation of rows and columns. We use 
the notation a,b, . . . for indices a,b, . . . € / and d, b ... for a,b . . . € /. The solutions (|2.4p 
provide an evaluation homomorphism of the infinite-dimensional Yangian algebra into some 
finite dimensional algebra. The latter we denote by 



Aj = [(J) ® H^ft ■ 



(2.9) 



There are 2 n homomorphisms of such form corresponding to all possible choices of the set / 
and every representation of Ql(I) in Ai provides a distinguished representation of the Yangian 
algebra. When the set / is chosen to be / = {1, . . . , n} we find the well-known result 



= L il,..,n} W = * I + 4® Ja , 



(2.10) 



where by we denote fll(n) generators in the fundamental representation and J£ generate the 
jjt(n) algebra 

[J^J2 = 5 c b J a d -5 a d J c b . (2.11) 

While the Lax operators L/(z) are sufficient to construct Q-operators for the fundamental 
representation, the first step in the generalization of [2] is to find "R-operators for Q-operators". 
These operators are defined as intertwiners between a given representation of the algebra Ai 
and the chosen representation of the gl(n) algebra related to V . However, it is convenient not 
to specify these representations and work with the abstract generators of Ai in (|2,4p and gl(n) 
in (I2.10p . The defining relation for R-operators is then the YBE of the form 

C(z 1 )L I {z 2 )TZ I (z 2 - zx) = 1li(z 2 - zx^iiz^dizx) . (2.12) 

The intertwiners TZj(z) are the basic building blocks of the transfer matrices of the form 

Tr/ [X>j Ki{z) ®Ki{z)®...® Kr{z)\ . (2.13) 

The precise form of the regulator T>j and the definition of trace will be explained in the following. 



1 We are changing notation with respect to [2] as follows: 
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2.2 New Solutions of the Yang-Baxter Equation 



This section contains the derivation of the intertwiner TZi(z) entering (|2.12p . We will use some 
important results of classical invariant theory. The latter appear to be closely connected with the 
theory of Yangians, see |23| for more details and references to the large literature on the subject. 
In order to solve (|2.12p for TZj(z) it is convenient to make use of the difference property of R- 
operators. Using the explicit form of the solutions (|2.4[) and (|2.1U|) we obtain the commutation 
relations 

[4,n I ] = 4n I , [4,n I ]=n I 4, [n^jg + h%] = o. (2.14) 

Without loss of generality we may rewrite TZi(z) in the factorized form 

Tll{z) = e*" Jc c U ,i(z) e~ a i J " . (2.15) 

Remarkably, (|2.14|) immediately implies that Hoj(z) does not depend on oscillators and that 
it is Qi(I) invariant. In addition, using (I2.15p . we obtain n x n equations for TZqj(z) from the 
YBE, naturally organized into four blocks: 

[Kajizljf] =0, (2.16) 

K ,i{z)4 = [{z - |) J| + J h c Jt - 4 4) KoM > (2-18) 

Ho tJ (z) ({z - % ^ + J b d JS - 4 4) = ((* ~ f) Jb + 4JS ~ 4 4) KoM . (2.19) 

Let us focus on the case relevant for the construction of the Q-operators. This means that we 
restrict to, cf. (|2.9p . 

Ai = ■ 8) H {IJ) , (2.20) 

where we take the singlet representation, denoted by the dot •, of the subalgebra Qi(I) - we 
call (I2.20P a minimal representation. This means that the generators J in (|2.5p act trivially on 
(i.e. annihilate) all states of a given minimal representation. We will nevertheless keep equations 
(I2.16p - (|2.19p in their full form because we will need them to prove the fusion relations in the 
following section. In the case of minimal representations (|2.19p follows directly from (|2.17p and 
(|2.18p , and the solution of (|2.17p satisfies also (|2.18p . In the following we focus on (|2.17p under 
the restriction that TZqj(z) can be expressed as a function of the Ql(I) Casimir operator^]. 

As a warm-up exercise, and in order to investigate what types of solutions may be expected, 
we solve (|2,17p for the simple case when \I\ = 1. As the dotted indices take only one value we 
can simplify it to 

n 0J (z) {z + \-4) 4 = 4n ,i(z). (2.21) 

2 We believe that this condition follows from 

{■Roj , J?] = , [TZo,! , 4 ] = , [Wo,/ , 4 J% } = , 

which are specializations, respectively, of the last equation in (|2.14l) . (I2.16|) and Q2.19P in the case of minimal 
representations. Clearly IZoj can be considered as a function of the Casimir operators of gl(n) as well. These are 
just constants in a given irreducible representation and will not enter the discussion regarding the determination 
of Uaj- 
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In this case Hqi(z) can be expressed as a function of the linear Casimir J? of Ql(I). Now, using 

[4,41 = 4, => 4f(4) = n4-i)4, (2.22) 

we can turn ()2,2ip into a simple system of first-order recurrence relations. These relations may 
be solved up to a spectral parameter-dependent periodic function in J? of period 1. The solution 
of (|2,2ip can thus be written with the help of the Euler gamma function as 

%w = p(z, 4) n* + \- 4) • (2.23) 

Next, we observe that a similar factorized form of (|2.2ip holds also in case the generators of 
the algebra may be written in Schwinger form as Jg = 8l A 8lb or Jg = c A Cs, where, respectively, 
bosonic or fermionic harmonic oscillators are employed. This is always possible for representa- 
tions whose Young diagram is of one-row or one-column form. We remark that these types of 
representations are relevant for the integrable spin chain that appears in the spectral problem 
of planar Af = 4 Super Yang-Mills. One again finds that the solution to the YBE is expressed 
as a gamma function depending on central elements. We will see in the following that gamma 
functions also provide a solution for more complicated cases. 

After these initial considerations, let us next proceed to the generic case. Here I is any set 
with q = \I\. We again look for a solution of (|2.17p depending only on the Casimir operators of 
fll(I). Then both sides of (|2.17f> can be seen as elements of the space Xj spanned by 

(J k )t = 4 1 4l---4 k _ 1 , A; = 1,2,... , (2.24) 

with coefficients in the center of the universal enveloping algebra U (q 1(1)). An important ob- 
servation is that these coefficients have very non-trivial commutation relations with elements 
of Xj. It is therefore not straightforward to solve (|2.17p . In the above warm-up exercise the 
key simplifying feature was that the space Xj is one-dimensional. In the generic case, thanks 
to the non-commutative Cayley-Hamilton theorem reviewed in Appendix lAl the space Xj has 
dimension q - it is finite dimensional. This may be easily seen by multiplying the relation (|A.8D 
(with dotted indices in place of the capital ones) from the left by j\ 

(J Q+ X= E (j k Y a k (h,...,£ q ). (2.25) 

k=l ° 

The quantities £i,...,£ q play an important role in the following. The reader may find more 
details in Appendix lAl and Section \2. 41 In analogy to (|2.25p one can reduce any higher power of 
J and express all of them in the basis given by 

k = l,...,q. (2.26) 

However, it turns out that this is not a very convenient basis. We are rather looking for a basis 

(X k )\ k = l,...,q, (2.27) 

such that for any element a L in the center of U(gl(I)) there exists another element a R in the 
center of U(q\(I)) for which 

a L [x k Y = [x k Y a R . (2.28) 

Finding such a basis is equivalent to the diagonalization problem of a family of commuting qx q 
matrices, which simplifies to the diagonalization problem of only one such matrix. The simplest 
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non-trivial element of the center of U(q\(I)) to be considered is the quadratic Casimir operator 
C*2 = J\Jf- One easily derives its action on the basis (|2.26p 

C 2 = 2 (J k+l ) a a + (J k ) a a (C 2 + q) ■ (2.29) 

Thanks to (|2.25ft . the determination of the basis f|2.27[) reduces to the diagonalization of the 
single matrix 2 M + (C2 + q) I, where 



1 


.. 


. 


o,\ 


\ 


1 


.. 


. 


«2 







.. 


. 






V 


.. 


. 1 


a q 


J 



(2.30) 



The elements a\ in the last column of M are symmetric functions of the £i appearing in (|2.25p . 
Their explicit form is given in (|A.9|) . We immediately conclude that after diagonalization one 
has 

6 2 (* fc ) ° = ° (62 + q + 2 4) , (2.31) 

where the l\~ are elements of the center of U(qI(I)) (see formula (|A.4|) ). We stress that up to an 
overall normalization the basis X k is uniquely fixed once M is diagonalized. This means that 
the basis elements X k have good exchange properties with all higher Casimir operators, and 
therefore with any element of the center of W(g [(/)), and thus all £i. The explicit form of the 
generators in (|2.26|) in terms of the basis (|2.27|> is given by (lA.lOp . 

Let us elaborate on the role of the 1^ and their relations to the £)[(/) Casimir operators 

C i = J£J%...J%_ 1 . (2.32) 

The operators t{ act by scalar multiplication on the elements of any $1(1) irreducible represen- 
tation. For highest weight representations their action may be determined by acting with the 
g[(/) generators on Ql(I) highest weight states 

JJ \h.w.s.) =0, for 6 < d , and J| \h.w.s.) = Ac \h.w.s.) . (2.33) 

In this case the operators 1^ act as shifted weights, and we deduct li\h.w.s.) = (Aj — i+l)\h.w.s.}. 
On a given irreducible representation of gl(n) the Ql(I) shifted weights 1^ should be understood 
as operatorial shifted weights. Upon acting on a highest weight state as in (|2.33p any polynomial 
in the Ql(I) Casimir operators can be rewritten as polynomial in the shifted weights £j. This 
procedure provides an isomorphism between the center of Ql(I) and the universal enveloping 
algebra of the Cartan subalgebra of . This is known as the Harish-Chandra isomorphism, 
see for example |23| and references therein. In particular one gets (see |24| ) 

& = e n ( 1 + i~zt) {ik)l ■ (2 - 34) 

k=l j^k \ lk l o J 

It allows us to write down the commutation relations between X k and l{ operators as 

k (x k Y c = [X^ (5 ik + , (2.35) 

3 For convenience we rewrite \t i as Ai. This rewriting manifests the fact that the natural ordering in the set 
{1, 2 ... ,71} induces an ordering on the set I. The same ordering is used in (|2.33p . 
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which is consistent with (|2,3ip . The exchange relations ([2.35|) are entirely fixed considering the 
analog of (|2.3ip for higher Casimir operators together with (j2.34j) . 

We are now ready to solve ([2.17p . We use the relation f|2.35[) and the relation between J k 
and X k , which we established in the process of diagonalization of the matrix M. This leads to 
the set of equations 

K 0)I {z) [z-\- £{ + 1) (x k ) C c = [X^ K j(z) , (2.36) 

which are direct analogs of the |7| = 1 case (|2.2ip . We added an extra label / to £^ in order 
to emphasize that these are shifted weights of the subalgebra qI(I). The solution for 1Zqj(z) is 
then given by 

<? 

n ,!(z) = Pl (z) [] F(z - | - £{ + 1) . (2.37) 

fc=l 

The quantity p(z) is an arbitrary function of the spectral parameter. Moreover, it may be an 
arbitrary periodic function of £^ of period one. We stress that (|2.37p is a symmetric function 
in the variables £[,... ,£ q . This is no coincidence since the Weyl group of gl(I) acts on these 
variables by permutations. The expression for 1Zqj{z) is universal and independent of the 
existence of highest (or lowest) weight states. 

To summarize, ()2.37p together with the factorized expression (|2.15p yields the solution to 
the YBE (|2.12p . We thus found the building blocks for the Q-operators for the spin chain with 
any Q\{n) representation A in quantum space. Before proceeding to the explicit construction of 
the corresponding novel Q-operators we will first prove some very important fusion relations. 

2.3 Fusion of Lax Operators 

In this section we present the factorization formula of the Lax operators that leads to the 
full hierarchy of functional relations. This procedure was exhaustively discussed in [2j for the 
fundamental representation of Ql(n) in quantum space. It provides a solid starting point for 
the factorization, based on a precise relation between certain degenerate representations of the 
Yangian algebra. In particular, in the present language, the structure of the auxiliary space is 
independent of the one of quantum space. Therefore the same functional relations remain valid 
for arbitrary representations A. Many more details and results on this issue are provided in |25| . 
In generalization of [2] we obtain the fusion relation for non-intersecting sets I T\ J = as 

Hf{z + A + f ) Tlf(z - f ) = Snfjj(z) 4 2 '1 S-\ (2.38) 

with 

4 2 ' ] =ex P (-^(a[ 2 '])'), (2.39) 

and S given by formulas (3.8)-(3.10) in [2]. In this section dotted indices refer to the subset J, 
and we define p\ = \I\, p2 = \J\ and p = pi + P2- The similarity transformation S acts purely in 
auxiliary space and disentangles the oscillators, making manifest that IZnjj(z) and TZq act in 
different auxiliary spaces denoted by [1'] and [2'], respectively. After a careful analysis of (|2.38p 
one concludes that IZiuj(z) is of the form (|2,15p with 

<L = e^' ] ^^o,K^i)e-( a[1 ' I ^^o, J (. 2 )e-^' I ^ C ? . (2.40) 

Via tedious but straightforward calculations one sees that 1Zo,luj{ z ) satisfies (|2.17p with 
in the Holstein-Primakoff form built from oscillators a^ 1 1 and a^ 1 1 (see formula (3.12) in [2] for 
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their detailed form). It means that TZjuj(z) is a solution of the YBE (|2,12|) with the g[(/ U J) 
representation in Amj given by weights A = C^^^O^,)- 

PI P2 

Following the same line of reasoning as in |2], we can summarize this section with the 
factorization formula for Lax operators 

K I {z)K {h} (z + X 1 )K {i2} (z + >! 2 ) ...K {lp} (z + \' p ) = S I Kf(z\A p )K Gl Sy 1 , (2.41) 

where / = {ii, . . . ,i p } and Aj = Xj + p ~ 2 ^ +1 , and kj depends on the normalization of 1Zj{z). 
Here we introduced an infinite dimensional highest weight representation 7T^ of the algebra 0l(p) 
with weights A p = (Xj)j=i p and 

1l+(z\A p ) = Tr+ p [K I (z)}. (2.42) 

The form of (|2,4ip is exactly the same as the one obtained for the fundamental representation in 
quantum space. As in the latter case, when supplemented by the known relation between finite- 
dimensional and infinite-dimensional modules (the so-called Bernstein-Gelfand-Gelfand (BGG) 
resolution) it suffices to derive the full hierarchy of functional relations. This means that the 
functional relations will be exactly the same as in the earlier case. The difference appears merely 
in the analytic properties of the Q-operators, which are a direct consequence of formula (|2.37p . 
In the following section we shall take a closer look at the analytic properties of TZi(z). 

2.4 Analytic Properties of Lax Operators 

Before we proceed with the construction of Q-operators for general lowest /highest weight gl(n) 
spin chains we would like to investigate the analytic structure of the new solutions to the YBE 
as found in Section \2. 2 1 This analytic structure will be reflected in the form of the Q-operators, 
which will be crucial in Section 12.61 We want to stress that it follows directly from our opera- 
torial construction of Q-operators. Three types of £)l(n) representations A will be distinguished 
in quantum space: compact representations, non-compact representations with a lowest /highest 
weight state, and non-compact representations without a lowest/highest weight state. The dif- 
ference between these three cases is reflected in the structure of the eigenvalues of the I operators 
entering (|2.37p . Their precise spectrum crucially depends on the representation under study and 
can be extracted in a case-by-case analysis. However, a common feature of all representations 
is that the eigenvalues of v for a fixed set / are separated by integers. 

For the compact case the spectra of all t are bounded from both above and below. This means 
that there exists a normalization such that all of the TZi(z) have polynomial matrix elements. 
This information is essential when deriving Bethe equations. For non-compact representations 
with a lowest /highest weight state some of the £ for a fixed set / are not bounded from above. 
Because of (|2.37p . it is not possible to find a normalization for which the TZi(z) is polynomial - 
its matrix elements are meromorphic functions. Let us elaborate on this point with the simple 
example of the spin — \ representation of 0l(2), which is a unitary representation of u(l, 1). We 
may realize it with the use of two oscillators, and define a basis in representation space by 

\M) = (ab) M |0), M = 0,1,.... (2.43) 

The generators of fll(2) can be written with use of oscillators a, b aa3 

j{ = aa , J\ = -ab , (2.44) 
J\ = ab , Jf = -1 - bb . (2.45) 

4 Algebraically, this could also be rewritten by a "particle-hole transformation" in the standard form j\ — a 6 a a . 
To get the appropriate real form u(l, 1), however, we should use the above form, and impose the reality conditions 
a = and b = b^ . 
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^{1,2,3} 





^{1,3} 




X 



X 




[' 





■0 



Figure 1: Hasse diagram for the unitary highest weight representations of u(l,2) algebra. The 
brackets [ ] indicate that the enclosed Lax operator has non-polynomial matrix elements. If 
the bracket is absent, one may normalize the Lax operator such that its matrix elements are 
polynomial. 



We have to consider two cases: / = {1} or I = {2}. In both of them we are interested in 01(1) 
subalgebras of gl(2). We notice that fW|M) = J\\M) = M\M) is not bounded from above 
while l^\M) = jj\M) = -(1 + M)\M) is. This means that we can find a normalization for 
which 1Z{i}(z), depending on £^ 2 \ t is polynomial. This is not possible, however, for 7£{2}(<z), 
which is a function of ^W. 

Considering the appropriate reality conditions, a similar analysis can be applied to oscillator 
representations of u(n,m), or even more generally to any unitary highest weight representation 
of u(n,m). The case of u(l, 2) is given in Figure [TJ where all 7£-matrices are depicted on the 
Hasse diagram. Similar patterns can be obtained for a generic highest weight infinite-dimensional 
representation of gl(n). For representations without the lowest/highest weight state none of the 
t operators is bounded from above or below and therefore none of R-matrices can be polynomial. 



2.5 Construction of Q-operators 

Let us now proceed with the construction of the Q-operators for arbitrary representations of 
gl(ra). We closely follow [2j, where the trace-construction was studied for the fundamental 
representation of gl(n) in quantum space. There we had introduced more general chain operators 
X/(z), whose building blocks were the Lax matrices in (|2.4p . The notion of X-operators can 
be easily generalized to the case of an arbitrary representation in quantum space by the simple 
substitution L/ — > IZi. To keep formulas short, we exclusively focus in the following on the case 
of Q-operators. 

The quantum space V of the fll(n) spin chain with the local state space V being in the 
representation A is an L-fold tensor product space, cf. (|2.ip . To construct the Q-operators we 
consider the monodromy matrices 



M/(z) = Vi Hi(z) ® Hi{z) ® . . . ® Ki(z) 



(2.46) 



acting on 




(2.47) 
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The local building blocks of the M/ are the R-matrices TZi(z) given in (|2,15|) . which were 
obtained from the YBE (|2.12p . While the various R-matrices in (j2.46|) act on distinct spaces V 
they all share the same auxiliary space Ai- The boundary twist operator T>j acts only in the 
auxiliary space. Its explicit form is given by 

X>/ = exp{-i*>hg}, (2.48) 

where = <£ a — $^ and h? = a^a? + | (no summation). The Q-operators are defined by taking 
the traces of the monodromy matrices (j2.46j) as 

Qfe) = e te ^. €J *-)Tf w(JiJ) { Mj (2r)} , (2.49) 
where, following the notation of [2], the normalized trace is defined by 

Tr n {e~^ h * } = Tl ^r $ LV • ( 2 -50) 

An important feature of the normalized trace is that it is completely determined by the commu- 
tation relations (I2.6P and the cyclic property of the trace. Recall that, thanks to the underlying 
YBE, all Q-operators above form, for a fixed quantum space, a commuting family of operators. 
They belong to a larger family of commuting operators X(z), which includes the standard trans- 
fer matrices as special cases. The next task is to diagonalize the Q-operators. This is done using 
the QQ-relations, to be derived in the next section, in conjunction with the analytic structure 
of the eigenvalues of the Q-operators. 



2.6 Functional Relations and Bethe Equations 

Applying the reasoning from [2j to the fusion relations of Section 12.31 one obtains two impor- 
tant relations between Q-operators, allowing to diagonalize the problem with the use of Bethe 
equations. All Q-operators can be expressed using only partonic Q-operators. These correspond 
to the sets / = {i}, containing only a single index i = 1, ... ,n. The formula is given by the 
following determinant^ 

Ai(*)Qj(*) = det||Q ai (z-j + 2±i)|| 1 < iij <| / |, (2.52) 
where Aj depends only on the twist angles defined in Section [2.51 It reads 

A/($)= 11 2i sin . (2.53) 

l<i,j<|/| 

As a direct consequence of these relations one can write down the QQ-relations, see also [26] . as 

A {a ,6}($) QlUaUb(z) Ql(*) = QlUa(* ~ \) Qlub(z + \) ~ Qlub(z ~ 5) QlUa(z + |) . (2.54) 

Just as in the earlier case of the fundamental representation in quantum space we can associate 
all Q-operators to the vertices of a Hasse diagram, which for gl(n) takes the form of an n- 
dimensional hypercube. Each of its faces corresponds to precisely one relation (|2.54|) . With 

5 One can check that this holds if we take pi(z) — 1 in (|2.37|) for all I. Another interesting normalization is 

n / n \ _1 

pi(z) = n r ( z + 1 - n r t z - ii - 1/ 2 +i n ' ( 2 - 51 ) 

i=q+l \i=l ) 

in which case all Q-operators are polynomial for compact representations. 
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every path along the edges of the diagram leading from Q (z) to Qii !M .. >n \{z) we can associate a 
system of Bethe equations. The different paths correspond to the different sequences of ascending 
sets = Iq C I\ ■ ■ ■ C I n = {1, 2 . . . , n}. Clearly there are n\ such paths. The functional form 
of Bethe equations is then given by 

1 Qu-A*i-h) Q^+i) Qu+M'-h) , 2 55) 

where zj* are the roots of Qj^z), and we denoted the eigenvalues of the operator Qi(z) by 
Q/ !,--)• 

So far we did not use any information about the analytic structure of the Q-operators, apart 
from the fact that there exists at least one root of Qi(z). It is needed in order to rewrite (|2.55D 
into the standard explicit form of Bethe equations. This is possible whenever all operators Q/ 4 {z) 
for i = 0,...,n on the chosen path can be normalized to be polynomials (see the discussion in 
the Section \2A\) . Let us start from compact representations in quantum space. Then for all 
Q/(z) we can write their eigenvalues in the form 

(q \ L m/ 

a=l J k=l 

Substituting this form into (|2.55p we end up with 

/ Ji _ n-i-l _ f _ l \ L 

e *(*a i+ l-*ai) I 2 2 = 

\ Ji _ n-i-l _f I 1 I 



(2.57) 



n z f ' z fc' 2 TT Z l [ Z k + 1 TT fL 
y Ji /<-i , 1 11 _A _ _-f< _ i 11 J 



+2 k+l Z l Z k 1 k Z k Z ^2 

These equations appeared for the first time in jB], where they were derived by the algebraic 
Bethe ansatz technique. For compact representations, due to Weyl permutation symmetry, 
we have that for fixed |/| all Qj(.z) are isospectral. This means that all n! systems of Bethe 
equations have exactly the same form. This is not true for non-compact representations - Weyl 
symmetry is broken, a fact which leads to different analytic properties of the Q-operators. For 
representations with a lowest /highest weigh state, however, there still exist paths on the Hasse 
diagram for which all Qi(z) connected by the path can be normalized to be polynomial. In that 
case the formula (|2.56p holds true and we are again able to rewrite (|2,55p into the standard form 
(|2.57p . This means that, using the Q-operator approach, for lowest /highest weight compact or 
non-compact representations it is always possible to write polynomial Bethe equations for spin 
chains. 



2.7 Q-operators in use 

In the previous sections we obtained new solutions of the YBE and used them to construct 
Q-operators. We furthermore proved operatorial relations between them, and exploited these 
to derive Bethe equations without the use of a Bethe ansatz. However, for practical purposes 
such as finding Bethe roots or eigenvectors from Q-operators, it is necessary to find their exact 
operatorial form by constructing suitable monodromy matrices and explicitly evaluating the 
trace over their auxiliary space. With increasing length of the spin chain and dimension of the 
quantum space this becomes a very cumbersome problem. 
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For compact representations the Q-operators are (dim\) L x (dim\) L matrices, where diniA 
denotes the dimension of the representation A, and L is the length of the spin chain. We 
will present some powerful Mathematica™ code which allows to construct concrete matrix 
representations of Q-operators for various specific compact representations. It is based on the 
Gelfand-Tsetlin basis, and serves as a checking device for many properties presented earlier. It 
may also be used to extract Bethe roots and eigenvectors for given representations. This code 
can be found in Appendix [Bj 

For non-compact representations with a lowest/highest weight state the construction of Q- 
operators is much more involved. The reason is that the representation space is infinite dimen- 
sional and therefore the size of the matrices becomes infinite dimensional. It is nevertheless 
sometimes possible to obtain explicit results, and to evaluate the trace analytically. For exam- 
ple, for the L = 1 spin chain for any 0l(2) representation A in quantum space the Q-operators 
can be written as 

Q {a} (z) = T(z -4 + 1) 2*1 (4 -£i,J% z + 4: 1 _J^ a) ) > ( 2 -58) 

where {a} U {d} = {1,2}, and 2*1(0, 6; c; z) stands for the hypergeometric function. One can 
obtain the same result solving directly the Baxter equation (see e. g. |27|). For longer chains it is 
possible to write Q-operators purely in terms of flt(n) generators. However, we are not aware of 
any closed formula in that case - they are written as infinite sums. It is nevertheless possible to 
extract at least parts of the Q-operators by using the fact that they are block-diagonal matrices 
with finite-dimensional blocks. 

3 Hamiltonians 

3.1 Stating the Problem 

The purpose of this section is twofold. Firstly, we provide an explicit and computationally very 
efficient form of the nearest-neighbor Hamiltonians for those integrable spin chains where the 
local representation space V in (|2.ip is in a so-called generalized rectangular fll(n) representation. 
It generalizes the famous digamma formula for the n = 2 case |10| as derived from the 3 1(2) 
R-matrix of Kulish, Reshetikhin and Sklyanin [9], which is widely used in many areas of the- 
oretical physics such as the AdS/CFT correspondence, QCD high-energy scattering problems, 
and integrable lattice models of condensed matter theory. Secondly, we show that the method 
we used in Section [2] for finding the Lax operators for Q-operators is rather universal, and can 
be used also in other cases. A particularly interesting feature of our result is the appearance of 
operators corresponding to shifted weights of tensor product representations. 

In Section [2] we used the QISM as a guiding principle, and derived the solution of inte- 
grable nearest-neighbor spin chains for arbitrary lowest/highest weight representations of fll(n) 
without employing the Bethe ansatz. This results in (equivalent) sets of Bethe equations, and 
simultaneously diagonalizes a large class of operators: standard transfer matrices T(z), Baxter 
Q-operators Q(z), and mixed forms X(z). In [12J it was shown that integrable nearest-neighbor 
Hamiltonians belong to the same commuting operator family. As such they are a consequence 
of the YBE, and thus naturally embedded in the framework of the QISM. To be more precise, in 
order to extract the Hamiltonian one has to construct the transfer matrix operator Ta(z), where 
the representation of the auxiliary space A is identical to the one of the local one-site quantum 
space V of ([2.ip . This special transfer matrix contains the information about all conserved local 
charges such as the shift operator, the nearest-neighbor Hamiltonian, and a further tower of a 
total of L higher conserved local charges. In order to extract all of these, one has to expand 
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the transfer matrix in the spectral parameter z around a special regularity point z = z*. One 
then reads off the various local operators as the coefficients of a Taylor expansion around z*. 
The regularity point z* is defined by the requirement that the building block of Ta(z), namely 
TZaa(z), turns into the permutation operator at z = z* as 

K AA {z*) = P. (3.1) 

In general, the process of finding TZa\(z) reduces to solving the appropriate YBE. When the 
representation A is fixed it is sufficient to consider the YBE on the tensor product space □ giAl^A, 
where we denoted the fundamental representation of Ql(n) by □. For generic representations 
solving the YBE is a tedious task. However, one may use a shortcut, exploiting the fact that 
TZ\\(z) is a A ® A-invariant operator. This means that the following equations are satisfied 

[^AaO), J h a ®l + l® 4} =0, a,b = l,...,n, (3.2) 

where J% in both spaces are evaluated in the representation A. For compact representations, 
for instance, there is a finite and rather small number of independent numerical matrices Invj 
satisfying relation (|3.2p . Then the R-matrix can be written as 

where ai(z) are spectral parameter-dependent functions. Substituting this form of R-matrix 
into the YBE immediately fixes the functions ai{z). This way one can find TZ\\(z) and thus 
the Hamiltonian density in terms of a series for any fixed representation A in quantum space. 
However, the question remains whether there is some closed- form expression for TZ\\(z) which 
works for any representation, i.e. is written purely in terms of the abstract algebra generators. It 
is known in the literature that the answer is partially affirmative for the case of those representa- 
tions whose tensor product decomposition is multiplicity free [15] • In that case, using the "tensor 
product graph method", one can write TZ\\(z) as a sum of projectors, with the coefficients de- 
pending on the second Casimir of the representations in the tensor product decomposition of 
A (Ei A. However, the form of TZaa( z ) obtained by using this method is not explicit, as one needs 
to construct case-by-case the needed projectors. 

The goal of the present Section [3] is then to provide a completely new form of TZ\\(z) for all 
rectangular (and therefore multiplicity-free) representations. The new form is both practically 
efficient as well as suitable for abstract considerations. The resulting R-operator is written in 
terms of gamma functions of the shifted weights of the tensor product algebra. When expanded 
in the vicinity of the permutation point it gives a general form of the Hamiltonian density as 
the sum over digamma functions, in generalization of the gl(2) case. In order to find the general 
form of TZ\\(z) we proceed as in the previous section, and solve the YBE equation without ever 
specifying the representation A. The solution is, up to a trivial multiplicative factor, unique, 
and is in perfect agreement with the one found using the tensor product graph method. 

In order to connect back to Section [2j we point out that it is also, in principle, possible to 
obtain the transfer matrix T\(z) by means of suitable fusion of Q-operators. This furnishes 
an alternative in-principle method for deriving the nearest-neighbor Hamiltonian of any spin 
chain with given representation A of V in (|2,ip . While we do not know how to directly find a 
"nice" formula for TZ\\(z) from this procedure, fusion nevertheless allows us, using functional 
relations, to write the dispersion law relating the set of Bethe roots found from the system of 
Bethe equations to the eigenvalues of the Hamiltonian. The derivation of the dispersion relation 
along these lines is presented in Section 13.41 
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3.2 Yang-Baxter Equation, Revisited 

We focus in this section on the jj[(n) representations A satisfying the generalized rectangularity 
condition 

J c b J a c =aJZ + p5 a h . (3.4) 

This condition is always satisfied for the gt(2) generators as a consequence of the non-commutati- 
ve Cayley-Hamilton theorem. For n > 2 the condition (I3.4|) imposes restrictions on the class 
of representations considered. Such representations of fll(n) shares some important feature with 
representations of Ql(2), most notably the tensor product of two such representations as well as 
their weight diagram are multiplicity free. For finite dimensional representations one can check 
that (3.4) is satisfied if and only if the corresponding Young diagram is rectangular. In this case 
a = s — a and f3 = sa, where s is the number of columns and a is the number of rows in the 
Young diagram associated with A. In general one can use the Holstein-Primakoff representation^] 
to prove that (j3.4|) is true in the case 

Ai = A2 = • • • = A a 7^ A a+ i = ... = X n , (3.5) 

where A, are the Dynkin labels of the representation A (see also |28|). This is true for the 
rectangular compact representations where A a > A +i, and corresponds to non-compact repre- 
sentations in the converse case. 

As was pointed out in Section [3.1[ we would like to find the general form of TZ\\(z) when 
expressed in terms of the generators of gl(n). To this aim we solve the YBE 

C( Zl )C(z 2 )m(z 2 - z x ) = <R{z 2 - zx)C(z 2 )jO(zx), (3.6) 

where £{z) is given by (|2,10p . We have also replaced TZ\\(z) by 9\(z), since we would like to 
leave the representation A unspecified for the moment. Assuming P!SH(z)P = Jt(z), we can 
rewrite (13. 6D in the form 



<r h - r h ) + (« - J c b J a c)) £*(*) = k{z) [ Z {4 - 4) - {j c b r c - r b j a c )) , (3.7) 

where by J and J we denote the generators of the two copies of fll(n). We now solve this 
equation using a method very similar to the one in Section 12.21 the Cayley-Hamilton theorem 
again being the main ingredient. This time, however, we use it for the tensor product algebra 
instead of the subalgebra of fjl(n). We only sketch here the main steps of the derivation, which 
goes along the same lines as in Section 12.21 Let us introduce a basis by defining 

(3*+l)a := (JC! _ Jc 1)(J c 2 + j£ )(J3 + JC3 } (J a + jaj ( g >8) 

From the Cayley-Hamilton theorem we prove that the space spanned by the Z k is finite-dimensio- 
nal, and that 

n 

(r +1 ) a b =E (3 k Ya k (Ax,...,A q ), (3.9) 

k=l 

where Ai are shifted weights of the tensor product representation. After a straightforward 
calculation one finds the commutation relations of the basis operators (|3.8p with the second 
Casimir of the tensor product representation 



43 fc+1 -4a3 fc , (3.10) 



6 To be more precise these representations are exactly of the type entering (|2.40p in the case in which 1Zo,i(z) 
and TZo, j(z) corresponds to minimal representations. 
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where a is the constant in the relation (|3.4p . Following the same line of reasoning as before we 
diagonalize the second Casimir C 2 and obtain 

C 2 X k = X k (C 2 + 4A k -4a), k = l,...,n. (3.11) 

This relation is solved by 

Aj,Xk = Sj h Xk - $jk%k , k = l,...a. (3.12) 

We now have all ingredients necessary to solve (|3,7p . Let as rewrite it first in the basis (|3.8p 

((z - a)tf)i + (3 2 )£) !*(*) = ((* + a)tfT b - (3 2 )£) . (3.13) 

After converting to the basis and using the fact the constant in (j3.4|) is given in terms of the 
tensor product representation weights as 

^4, _i_ A- — \ 

a = — 5 , with ife = 2a- + (3.14) 

we find 

X, f , + + SK(*) = <H(z)£ fe (z - Ak ~^ + 1 ^j , fe = 1, . . . , a , (3.15) 

The solution is then derived to be 

= f[ r (z + ifc ~^ + 1 l r ^ - Ak ~f + 1 + 1^ . (3.16) 

The function p(z) can be any function of the spectral parameter and any periodic function of the 
A. In Appendix ICl we present the check that this formula indeed agrees with the one found from 
the tensor product graph method. In the following we will further demand that the R-matrix 
satisfies the unitarity condition 

<R{z)<X(-z) = 1 . (3.17) 
This fixes the normalization of %\{z) to be 

P (z)= ^nr^ + Afc + 2 t + 1 + 2o - fe ) r ( z - i4fc + ^ + 1 - 2a + A; + 1 )) ' < 3 - 18 ) 

unique up to any function satisfying (|3.17p . In order to apply this result to a specific representa- 
tion A, we merely need to change the notation back from 9\(z) to TZ\\(z). Since the expressions 
(|3.16p . (|3,18p only depend on the shifted weight operators A k and the spectral parameter z, they 
may be immediately applied to any A. 



3.3 Permutation and Hamiltonian Density from the R-matrix 

We next use equation ([3.16p to derive an explicit form of the permutation operator and the 
Hamiltonian density for the spin chain. The result will again be written in terms of the shifted 
weights of the tensor product representation. The operators are given, respectively, by the first 
and the second term of the expansion of 9K(z) in the spectral parameter 

9t(z) =P(l-zH + ...). (3.19) 
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Then the permutation operator is simply 9t(0), which from (|3.16p and by using identities for 
the gamma functions reads 

a i . . 

p = nt -1 ) ■ ( 3 - 2 °) 

k=l 

In order to extract H we take the logarithmic derivative of JH(z) and evaluate it at the permu- 
tation point z = z* = 0, yielding 

U = — ^-logJK(z) . (3.21) 

z=o 

After straightforward calculations, which are presented in Appendix [Dl we end up with the 
following form of the Hamiltonian density for the spin chain with the representation in the 
quantum space satisfying the generalized rectangularity condition (|3.4p : 



k=l 



(3.22) 



It is easy to see that it reproduces the well-known |10| . |13j form of the Hamiltonian density for 
su(2) spin chain with spin s, namely 

U s = 2V>(s + 1) - 2tp ~2 2 + 1 ^j = 2h ^ ~ 2/l W- ( 3 - 23 ) 

Let us pause and make an interesting observation. It is clear from (I3.22p that the value 
of the Hamiltonian density for a given multiplet on the right hand side of the decomposition 
depends only on the shape of the Young diagram but not on the algebra at hand - it is the 
same for gl(2) and for g[(102). This means that the energies appearing in the spectrum of the 
periodic (in particular, all twist angles are zero) spin chain with length L = 2 are exactly the 
same for all gl(n). Dependence on the rank appears, however, when we calculate the dimensions 
of multiplets that grow with the rank of the algebra. This fact is not anymore true if we consider 
longer chains because in that case it is known that the spectrum cannot be fixed completely by 
symmetry arguments alone. 

We recall that, as in [2J, the nearest-neighbor Hamiltonian commuting with the family of 
transfer matrices constructed in the previous section is 

L 

H = J^ >i+ i, with Ww^e^^W^^^W, (3.24) 
i=l 

where "Hi,i+1 1S the Hamiltonian density given in (|3.2ip . The presence of the twist angles (f) a 
breaks the gl(n) invariance of the chain and allows to probe the "fine structure" of the spectrum. 



3.4 Energy Formula from Functional Relations 

The process of diagonalization of the Hamiltonian for the spin chain using Bethe equations is 
not complete unless we supplement them with the dispersion relation, which allows to determine 
the spectrum of the Hamiltonian in terms of the Bethe roots. We show here that for the finite- 
dimensional representations the operator Ta(z) generating the Hamiltonian can be expressed 
solely as a function of a chain of Q-operators connected by a single path in the Hasse diagram. 
Then, the eigenvalues of the former (energies) can be written in terms of eigenvalues of the 
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latter (functions of the Bethe roots). In order to extract energies from Ta(-z) we write for every 
eigenvalue 

E = -—logT A (z)\ z=Zi> . (3.25) 

We recall that the transfer matrix belongs to a larger family of operators, namely X-operators, 
and that for the latter we can write, again following [2j, 

A Ik Xj k (zi,...,z k )= det Q a {z b ). (3.26) 

In (|3.26p we omitted the explicit dependence on the spectral parameter, which may be recovered 
for every X/ fc (z) by taking the mean value of all the numbers z present in the expression. From 
the properties of determinants we can write 

Xj fc (z ,z 2 , ...,z k ) _ X Ik (z ,z 1 ,. . ■ X Ik _ 1 (z 2 , ■■■ ,Zk) X Ik _ 1 (z ,z 2 , ■ ■ • ,-Zfc-i) 

X Ik (z 1 ,z 2 ,...,z k ) X Ik (zi,z 2 ,...,Zk) X Ik _ 1 (z 1 ,z 2 ,...,z k _ 1) X Ik _ 1 (zi,z 2 ,...,z k _ l )' 

(3.27) 

which is is just a version of Pliicker's relations, and thus valid for any numbers z a . 

By an appropriate choice of the numbers z a we can produce any transfer matrix on the 
left hand side of (|3.27|) . Then the right hand side of (|3.27p will contain operators which are 
simpler in terms of the representation involved. Repeating this procedure, we can reduce all 
representations to the minimal ones. The right hand side will then only depend on the operators 
of the form Xj k (zo, ... ,zi) with z% — z%-\ = 1, which are exactly the Q-functions. Interestingly, 
we conclude that any transfer matrix is a sum of d terms, where d is the dimension of the 
representation of the auxiliary space of the transfer matrix. 

The analytic properties of all operators present on the right hand side of (|3.26p are fixed by 
(I2.56p . Furthermore, in line with our definition of Q-operators (|2.49|) . the regular point of Ta(z) 
is determined to be[j z+ = The procedure presented above is rather complicated for generic 
representations, and we have not followed it through in detail. It is however straightforward 
to check by studying small representations that the energy formula is as expected given by the 
expressio 




E -Z^2^ [ Ji n-i-l p — ~ J t n-^—p —J ■ ( 3 - 28 ) 

2 <-n—i+l 2 *l 2 ^n-i T 2 / 



A general proof is under construction. 
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A Review of Non-commutative Cayley-Hamilton Theorem 

We present here a non-commutative version of the Cayley-Hamilton theorem following [23j (chap- 
ter 7). It is used in Section f2.2l when solving the YBE for building blocks of Q-operators. It is 
remarkable that this theorem is naturally embedded in the context of Yangians. 

Let us first review the well-known Cayley-Hamilton theorem for matrices. Given a q x q 
matrix M with commuting entries one has 

p(M) = , where p(z) = det(z — M) = z q + a q z q ~ l -\ V a u (A.l) 

where p{z) is called the characteristic polynomial. The coefficients aj are symmetric functions 
of the eigenvalues of M. The statement above can be equivalently rephrased saying that the 
matrix M q can be written as a linear combination of matrices M k with k < q and coefficients 

Oj. 

Imagine now that we have a matrix whose operatorial entries satisfy the Ql(q) algebra. Let 
us define the non-commutative analog of the characteristic polynomial known as a Capelli de- 
terminant 

C(z) = £ (-If (* + ...(* + J - g + 1)^ • (A.2) 

We see that C(z) is by construction a polynomial in z with coefficients in the universal enveloping 
algebra U(gi(q)) and can be written as 

C(z) = z q + C q z q - 1 + •••+£!. (A.3) 

It can be proven that the coefficients C{ belong to the center of U(gl(q)). Moreover 

C(z) = (z + £ 1 ) ■■■ (z + £ q ), (A.4) 

where the Weyl group acts on l{ by permutations. The li are, up to a global shift, the same 
as the shifted weights Again, Cj are symmetric functions of the £. The non-commutative 
Cayley-Hamilton theorem states that 

C{-J T + q-l) =0, and C(— J) = , (A.5) 

where T means transposition as explained in the following. Powers of J should be understood 
as, e.g. 

{J 2 ) A B =J%4, A,B,C = l,...,q. (A.6) 

This means contraction of indices from top-left to bottom-right and analogously for higher 
powers. On the other hand, by powers of J T we denote the other contraction 

((J 2 ) T )i= JcJb, A,B,C = l,...,q. (A.7) 

Throughout this paper we use only the second formula of ()A.5P together with (1A.6|) . As for the 
commutative case, (|A.5j) can be rewritten as 

q-l . 

(J q ) A B=Yl (j k ) B a k+ i(h,...J q ). (A.8) 

fc=0 

9 For the convenience of the reader we point out that we will use a different looking, but equivalent use of 
transposition compared to [23] . This is done because we will use the second formula in (|A.5|) 



20 



We stress once more that a^i^i, ■ ■ ■ ,i q ) are symmetric functions of their arguments which can 
be written as 

a k (£ 1 ,...,£ q )= (-iy+ k Yl ( A - 9 ) 

ja<ji<-<j q -k 

We present here another formula we need in Section 12.21 

(JYa = A « W • A * W = 11^ - • (A.10) 

fc=l z^fc 
B Mathematica™ Code for Q-operators 

In this appendbJ^l we present a Mathematica code which allows to calculate Q-operators for 
compact representations. One has to specify the set I = . . . , i p } and the gl(n) representation 
A in the quantum space. The latter is of the form A={si, . . . ,s n } where Sk is the number of 
boxes in the k-th row of the Young diagram. The code is grouped in three partj^l: the first part 
calculates generators J* 2 of the algebra in the representation A with the function J J [i\ ,%i , A] . 
We also define Casimir operators a of the Ql(I) algebra with CC [i , I , A] . The second part defines 
the function Q[L,z,I,A] , which calculates the Q-operator Qi(z) for the spin chain of length L 
and representation A in the quantum space. Here z is the spectral parameter. In the third part 
we present some functions suitable for checking whether the Q-operators we constructed satisfy 
various basic properties such as proper commutation relations (using comm[L, setl , set2 ,rep] ) 
and the QQ-relations (using QQ [L, set , il , i2,rep] ). 



B.l Generators of Algebras 

ps [rep_] : =ps [rep] =Block [{res= (Prepend [rep , $0 &/@ (Tuples [Table [Range [rep [ [ 
l,i+l]] ,rep[[l,i]]] ,{i.l, Length [rep [[1]]]-1}]] ))}, If [Length [res [[1,1]]] ==1, 
Reverse@Sort@res , Reverse@Sort@Flatten [ps/ Ores , 1] ] ] ; 
dim[rep_] : =dim[rep] =Length[ps [{rep}] ] ; 

J0[k_] [pat_] :=J0[k] [pat] =Sum[pat [ [k, i] ] , {i , 1 ,k}] -Sum [pat [ [k-1 , i] ] ,{i,l,k-l}] ; 
Jp[k_] [pat_] :=Jp[k] [pat] =- Sum [Product [pat [ [k, i] ] -i-pat [ [k+1 , j] ] +j ,{j , 
Range [k+1] }] /Product [pat [ [k, i] ] -i-pat [ [k, j] ] + j , {j .Complement [Range [k] , {i}] }] 
i;[ReplacePart [pat , {k, i}->pat [ [k, i] ] +1] ] , {i , 1 ,k}] ; 

Jm[k_] [pat_] :=Jm[k] [pat] =Sum [Product [pat [ [k, i] ] -i-pat [ [k-1 , j] ] +j ,{j , 

Range [k-1] }] /Product [pat [[k, i] ] -i-pat [ [k, j] ] + j ,{j .Complement [Range [k] ,{i}] }] 

v [ReplacePart [pat , {k , i}->pat [ [k , i] ] - 1] ] , {i , 1 , k}] ; 

J J [k_ , 1_ , rep_] : = J J [k , 1 , rep] =Which [k==l , DiagonalMatrix [JO [k] /@ps [{rep}] ] , 
l==k+l, Transpose [Table [Coef f icient [Jp [k] [ps [{rep}] [[i]]] ,u[#]]&:/@ 
ps [{rep}] , {i , 1 , Length [ps [{rep}] ] }] ] , l==k-l , Transpose [Table [Coefficient [ 
Jm[l] [ps[{rep}] [ [i] ]] ,«[#]] &/@ps [{rep}] , {i , 1 .Length [ps [{rep}] ]}]] , l>k+l , 
JJ[k,k+l,rep] . JJ [k+1 , 1 ,rep] -JJ[k+l,l,rep] . JJ [k,k+l ,rep] ,l<k-l, JJ[k,k-l,rep] . 
JJ[k-l,l,rep]-JJ[k-l,l,rep] . JJ [k,k-l ,rep] ] ; 

CC [k_ , set_ , rep_] : =CC [k , set , rep] =Sum [Dot [SequenceSOTable [ J J [ j [i+1] , j [i] , rep] , 
{i,l,k}]]/. j [k+l]->j [1] , Evaluate [Sequence@@Array [{j [#] ,set}&,k]]] ; 
CCDiag [k_ , set_ , rep_] : =CCDiag [k , set , rep] =Diagonal [Inverse [Transpose [ 
Eigenvectors [CC [Length [set] , set , rep] ] ] ] . CC [k , set , rep] . Transpose [Eigenvectors [ 
CC [Length [set] , set , rep] ] ] ] 

10 TL would like to thank Nils Kanning for inspiring discussions about the code. 

n In order to use the code it is sufficient to copy it, and to paste it into a Mathematica™ notebook. 
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An example: the second Casimir of the gl(2) subalgebra of gl(4) generated by the set {1,3}. 
The representation of nt(4) is taken to be the fundamental representation rep={l, 0, 0, 0}. 
CC [2, {1,3}, {1,0, 0,0}] 

B.2 Q-operators 

Casimir0fL[p_,k_] : =Sum [Product [(1+1/(1 [i] -1 [j] ) ) ,{j .Complement [Range [p] ,{i}]}] 
Power [l[i] ,k] ,{i,l,p}] ; 

CHS [z_ , set_ ,rep_] : =Block [{n=Length [rep] ,p=Length [set] } , If [p>0 , (Transpose [ 
Eigenvectors [CC [p, set , rep] ] ] .DiagonalMatrix [Product [(Gamma [z-1 [il] - (p) /2+1] 
/. (Table [(NSolve [Table [CasimirOfL [p, i2] ==CCDiag[i2, set , rep] [[i4]] ,{i2,l,p}] , 
Table [l[i3] ,{i3,l,p}]] [ [1] ]), {i4, 1 .Length [CC [1 , set , rep] ]}] //Rationalize) ) , 
{il , 1 ,p}] /Product [Gamma [z- (rep [ [i5] ] -i5+l) -p/2+1] , {i5 , 1 ,n}] Product [ 
Gamma [z -rep [ [n] ] - (-i6+l) -p/2+1] , {i6 ,p+l ,n}] / /FunctionExpand] . Inverse [Transpose [ 
Eigenvectors [CC [p, set , rep] ] ] ] //Expand) .Product [Gamma [z-rep [ [n] ]-(-il+l) + l]/ 
Gamma [z- (rep [[il]] -il+l)+l] , {il , 1 ,n}] IdentityMatrix [dim [rep] ] //FunctionExpand] ] ; 
expp [set_ ,rep_] [A_] : =Block [{n=Length [rep] ,m=rep[[l]] , t},t=Flatten [Tuples [Tuples [ 
{Complement [Range[n] ,set] ,set}] ,#] &/@Range [1 ,m] ,1] ; Total [KroneckerProduct@@@({ 
Dot@0 ( J J [Sequence@@# , rep] &/@#) , Power [- 1 , Length [#] ] /Factorial [Length [#] ] D [A , 
Sequence@@#] &@ (x [Sequence@@Reverse [#] ] &/@#) }&/@t) ] +KroneckerProduct [ 
IdentityMatrix [dim [rep]] ,A]] ; 

expl [set_ ,rep_] : =Block [{n=Length [rep] ,m=rep[[l]] ,t} ,t=Tuples [{set .Complement [ 
Range [n] ,set]}] ; Dot@@( (IdentityMatrix [dim [rep] ]+Sum [Power [x [Sequence@@#] ,i]/i! 
Dot [Sequence@@ConstantArray [JJ [Sequence®^ ,rep] ,i]] , {i , 1 ,m}] )&/@t)] ; 
Lax [z_ , set_ , rep_] [A_] : =Block [{n=Length [rep] , p=Length [set] } , If [0<p<n , ( ( 
KroneckerProduct [expl [set , rep] . (CHS [z , Complement [Range [n] , set] , rep] / /N// 
ExpandAll/ /Chop/ /Rationalize/ /FunctionExpand) , IdentityMatrix [Length [A] ] ] // 
ArrayFlatten) . (expp [set, rep] [A] )//FunctionExpand//Simplify//Expand) ,CHS[z, 
Complement [Range [n] , set] ,rep] ] ] ; 

Q [L_ ,z_ , set_ ,rep_] : =Block [{a,p=Length[set] ,n=Length [rep] } , If [0<p<n, a=(Nest [ 
ArrayFlatten [Lax [z, set, rep] [#]]&, {{Product [Power [x [il , i2] ,nn[il,i2]] ,{il,set}, 
{i2, Complement [Range [Length [rep] ] , set] }]}} ,L] /Product [Power [x [il , i2] ,nn[il,i2]] , 
{il , set} , {i2 , Complement [Range [Length [rep] ] , set] }] / /Expand) / .Power [x [il_ , i2_] , i_ . 
] : >0//Expand, a=If [L>1 .KroneckerProduct [SequenceOOConstant Array [Lax [z , set , rep 
] [ {1}] ,L]] , Lax [z, set, rep] [{l}]]] ;Exp[I (z-Last [rep] )Sum[0[i] , {i , set}] ] a/ . { 
Power [nn[j_,i_] ,k_.] :> -(-l+Exp[I (0[i] -0[j] )] ) HurwitzLerchPhi [ 
Exp[I (<£[i]-<£[j])],-k,0]}] 

An example: Q-operator Q{x}(z) for the length L = 2 spin chain with the spin-1 represen- 
tation (rep={2, 0}) in the quantum space. 
Q [2, z,{l}, {2,0}] //Simplify 

B.3 Checks 

comm[L_,setl_,set2_,rep_] :=Q [L,zl , setl ,rep] . Q [L,z2, set2,rep] -Q [L,z2, set2 ,rep] . 
Q [L,zl , setl , rep] //Simplify; 

QQ[L_,set_,il_,i2_,rep_] :=Q[L,z+l/2,Union[set,{il}] ,rep] .Q[L,z-l/2,Union[set, 
{i2}] , rep] -Q[L,z- 1/2, Union [set, {il}] ,rep] .Q[L,z+l/2, Union [set, {i2}] ,rep]-2 i 
Sin [1/2 (</>[il] -0[i2] )]Q[L,z, Union [set, {il},{i2}] ,rep] . Q [L,z , set , rep] //Simplify ; 

An example: Commutation relation between Qi(zi) and (£2(22) for the length L = 2 spin 
chain in the adjoint representation of fll(3) (rep={2, 1,0}). 
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comm[2,{l} ) {2},{2,l,0}] 

An example: QQ-relation for the antifundamental g[(3) representation rep={l,l,0}. We 
take the length L = 2 and the plaquette containing operators Q{i}(z), Q{i,2}(- z )j Q{i,3}( z ) an d 
Q{l,2,3}(z)- 

QQ[2,{1},2,3,{1 J 1,0}] 



C Comparing R-matrix and Tensor Product Graph Method 

We want to prove that the formula (|3.16p which we found directly from the YBE reproduces the 
one found using the tensor product graph method (see e.g. |15|). Let us take any rectangular 
representation space V. Then the decomposition of the tensor product of V is multiplicity-free 
and reads 

v®v = ®v x . (C.l) 

We introduce a relation ~ for representations appearing in the decomposition (jC.ip as 

V\ ~ Vfj, V\ C V a dj <8> and V\ has different parity than . (C2) 

The relation ~ is symmetric and one can show that 

n 

v x ~v»^Y,\ A i- A i\ = 2 - (°- 3 ) 

1=1 

In words, it means that two representations are in the relation ~ if the former can be obtained 
from the latter by moving a single box in the Young tableaux. 
Prom the tensor product graph method we have 

Raa(z) = ]Tpa(^Pa, (C.4) 

A 

with 

pM= 6 —^, 6 = \(C 2 (\)-CM) ifonlyA~/x. (C.5) 
p^z) d + z 

Let us calculate the second Casimir 

ft - En (i - t^t) M- E (M - (* - + (2 ° - 2)( 6 2 ° " 1)2 ° ■ «*> 

i=l jyti V 1 3 / j=l 

Let us take A, fi such that A ~ fj,. Then using (|C.3P we have 

C 2 (A) - C 2 ( M ) = A\ + A\ - (i fc + l) 2 - (if - l) 2 = -2A k + 2A- k -2 (C.7) 

and 



Pa(^) = %(A- k -A k -l)-z 
P»{z) \{A k - A k - I) + z 
On the other hand from our formula we have 



r x ( z ) T { z + H A k -h + i)j r [z + |(A S - A fc + 1 

r (z + \{A k - A k + 1) + 1) r (z + \{A k - A k + 1) - 1 

z + \{A k -A k + l) _ \ (A k - A k + 1) - z 

z + UA- k - A k -\) UA k -A k -i) + z' 



(C.8) 

(C.9) 
(CIO) 
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The overall minus sign comes from the fact that in the tensor product graph method we consid- 
ered R(z) = PR(z). From the definition of the relation ~ we know that representations A and 
jj, have different parities then we can conclude that 

5^ = ^44, (en) 

proving the equivalence of the results (|3.16p and (|C,4p , 



D Derivation of Hamiltonian from R-matrix 



We use the formula (|3.2ip together with (|3.16p and the normalization (|3,18p in order to find an 
expression for the Hamiltonian density 



U 



d 

dz 



log 9*0) 



k=l 



1p\Z + 



2=0 

A k -A k + l 



A- k - A k + 1 



(D.l) 
(D.2) 



-ip\z + 



A k + A- k + 1 



+ 2a — k \ - ip \ z 



A k + A- k + 1 



2a + k + 1 



z=0 



One has to proceed carefully because some of terms in (|D.2|> are divergent while evaluated at 
z = 0. Still the result is finite because of the relation 



lim [ip(e — x) — ip(e — y)] = i]){x + 1) — ip(y + 1) , for x, y € N 

e— >0 



Using the fact that for eigenvalues of A we have 



Al - A k + 1 



6 



A k + A k + 1 



2a + k + 1 G 



(D.3) 



(D.4) 



we can write 



u 




4 



A k + A k + 1 



+ 2a-k 



ip(s + a- k + 1) 



(D.5) 
(D.6) 



where the last relation holds true for the compact representations with a rows and s columns. 
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